Decoherence-free generation of many-particle entanglement by adiabatic ground-state 

transitions 
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We discuss a robust and decoherence insensitive method to create many-particle entanglement 
in a spin system with controllable collective interactions of the Lipkin-Meshkov-Glick type and 
propose an implementation in an ion trap. An adiabatic change of parameters allows a transfer 
from separable to a large variety of entangled eigenstates. We shown that under certain conditions 
the Hamiltonian possesses a supersymmetry permitting an explicit construction of the ground state 
at all times of the adiabatic transfer. Of particular interest is a transition in a non-degenerate ground 
state with a finite energy gap since here the influence of collective as well as individual decoherence 
mechanisms is substantially reduced. A lower bound for the energy gap is given. 
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Entanglement is one of the most characteristic features 
of quantum systems and lies at the heart of quantum in- 
formation processing and computing 0. While in few- 
particle systems entanglement is by now reasonably well 
understood and practical schemes for its creation and 
manipulation are developed 0, Q , many-particle entan- 
glement is still an open field of research with a large unex- 
plored potential for applications. There exist for example 
proposals to implement quantum computation in an ini- 
tially entangled many-particle system by performing only 
measurements and single qubit operations, both of which 
are relatively easy to implement _.5;j. A necessary pre- 
requisite for this is however a specific many-particle en- 
tanglement. Since entangled states become increasingly 
susceptible to environmental interactions if the number 
of particles increases, an important practical challenge is 
the design of robust and most importantly decoherence- 
resistant mechanisms for its generation. We here pro- 
pose and analyse such a mechanism which is based on 
adiabatic transitions in a spin system with controllable 
collective interactions. Robustness against decoherence 
and short process times are achieved by choosing ground- 
state transitions with a large energy gap very similar to 
the ideas used in adiabatic quantum computation [(|. 

Let us consider a collection of N interacting spin 1/2 
systems described by a generalization of the Lipkin- 
Meshkov-Glick (LMG) Hamiltonian 
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where the Ji, (i € {x, y, z}) are the total-spin operators 
of the ensemble. We show that this Hamiltonian can 
be implemented for (i = with controllable parameters 
by a generalization of the ion-trap scheme suggested by 
S0rensen and M0lmer ||, where cold ions interact via a 
common trap oscillation and are driven by bichromatic 
lasers fields. In contrast to the effective Hamiltonian of 
the generalized LMG interaction can not be solved 
exactly. It does provide however the possibility for adia- 
batic and therefore robust transitions between separable 



and entangled many-particle states. The most important 
feature of is that such transitions are possible while 
staying in a non- degenerate ground state with a finite en- 
ergy gap and that a large class of entangled target states 
are accessible by varying fi. We show that in the special 
case of A = 1 , the LMG Hamiltonian has a supersym- 
metry and the ground-state can be explicitly constructed 
for all times of the transfer process. The gap can be made 
rather large and hence the process can be fast despite of 
its adiabatic nature. For the same reason the influence of 
decoherence is strongly reduced: Collective decoherence 
due to noise in the external control parameters is elim- 
inated by the adiabatic nature of the process and the 
influence of independent individual reservoir couplings is 
suppressed due to the presence of a finite energy gap. 



Let us first discuss the case of a negative coupling 
parameter £ < 0. If \i = and A > N, the system 
described by eq.Oj undergoes a quantum-phase transi- 
tion [j| when the interaction parameters \i an d X2 are 
changed (case i). This transition can be described analyt- 
ically in the semiclassical limit with J = N/2 and N 3> 1. 
If xi — X2, Jz is a conserved quantity and has a triv- 

ial anharmonic spectrum: H — > £xi [A J z — Jj + J ] . 
The ground state of the system has maximum total angu- 
lar momentum J = N/2 and a z projection m z = ±N/2 
depending on the sign of A. Both of these states de- 
noted by I ft . . . ) and | JJ . . . ) arc separable many- 
particle states. On the other hand in the limit Xi = 
the terms in containing J x and J z vanish and the 
Hamiltonian approaches H — > It has again a 

trivial spectrum, however with two degenerate ground 
states \m y = ±N/2). Both of these, taken individually, 
are separable. Symmetric or anti-symmetric superposi- 
tions of them form however the TV-spin analogues of the 
Greenberger-Horne-Zeilinger states (GHZ states). U sing 
the Schwinger-representation of angular momenta [lfj 
one finds that the state \m z = N/2) is adiabatically con- 
nected only to one particular superposition due to the 
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symmetry of the Lipkin interaction. Thus 
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which corresponds to the generation of the A-particle 
analog of the GHZ state in the a y basis. |±) denote 
single-particle eigenstates of a y with m y — ±1/2 respec- 
tively. Due to the phase factor e mJ the entangled state 
depends sensitively on the total number of atoms even in 
the limit N -» oo [Tlj. 

For a positive coupling parameter, £ > 0, three cases 
need to be distinguished: /i = and the total number 
of spins is odd (case ii); fi = and the total number of 
spins is even (case iif); and /i ^ (case iv). In all cases 
the ratio X1/X2 is again rotated from unity to zero. 

Since £ > the initial ground state in cases {if) and 
{Hi) is \m z — ±A/2) depending on the sign of A ^ and 
is separable. The final Hamiltonian is again H — £xl Jy, 
whose ground state is now however the eigenstate of J y 
with smallest value of \m y \. Thus for an odd number of 
particles (case if) there are two degenerate ground states 
with m y = ±1/2, both of them being maximally entan- 
gled. Making use of the symmetry of the interaction and 
using the Schwinger-representation we find that the adi- 
abatic transition leads to the mapping 
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where n = \N/2] and the subscript "s" denotes sym- 
metrization. 

In case {Hi) the final ground state (xi — ► 0) is non- 
degenerate, has spin projection m y = and is maximally 
entangled. 
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\m y = 0) = \{+) N / 2 {-m. (4) 



In this case there is no merging of eigenstates and conse- 
quently no phase transition. The absence of degeneracy 
during the entire adiabatic transfer makes case {iif) par- 
ticularly interesting because here decoherence is strongly 
suppressed by the presence of a finite energy gap. The 
target state is however fixed to the special case \m v = 0). 

The variety of accessible target states in a non- 
degenerate adiabatic ground-state transition can be sub- 
stantially increased by adding a linear interaction propor- 
tional to J y to the LMG Hamiltonian, i.e. by allowing 
for a non- vanishing value of /i in eq. Q , (case iv) : If we 
assume for simplicity that /i = to, with to being an in- 
teger or half- integer with m 6 {— J, — J + 1, . . . , J}, the 

final {xi — > 0) Hamiltonian approaches H — > 
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Its ground state is \m y — to) and thus can be ad- 



justed to have any eigenvalue of J y , which is maximally 
entangled unless to — ± J. 



In all four cases considered an adiabatic change of 
X1/X2 from unity to zero leads to the generation of a 
well defined entangled state. In cases {iif) and {iv) this is 
moreover possible with a finite energy gap. The adiabatic 
process is neither of the Landau-Zener nor the STIRAP 
type We have illustrated the different scenarios for 
even N, fi — 0, and a positive linear term in Q in Fig. 
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FIG. 1: Energy spectrum of Hamiltonian (0 for £ > 0, 
even number of particles N and maximum J in the limits 
for X1/X2 — > 1 and A > (left) and xi/xa (right). If the 
sign of £ is changed the picture has to be flipped upside down 
and if the number of spins is odd (case {ii)) the state m y = 
does not exist and the state \m y \ — 1 needs to be replaced 
by \m y \ = 1/2. Adiabatic change of X1/X2 allows for robust 
transfer between separable (left) and entangled eigenstates 
(right) following the possible scenarios {i-iii). 

To find the ground-state of for arbitrary values 
of X1/X2 as e.g. during the transfer process is a diffi- 
cult task. We will show now that possesses a super- 
symmetry (SUSY) for A = 1, allowing for an explicit 
construction of the ground state for cases {iif) and {iv). 
(The existence of an extra symmetry has been suspected 
[13j but has not been understood so far.) For A = 1 the 
Hamiltonian can be factorized as 

H 

T 



(xiJx + iX2Jy - «X2m) {xiJx - iXlJy + «X2^) 

(6) 
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Due to the SUSY the spectrum consists of twofold de- 
generate and a nondegenerate state 0|- F° r = rn £ 
{— J, . . . , J} the ground state is nondegenerate and obeys 
{XiJx - iX2J y + iX2m) \ip ) = 0. One finds 



\ip ) = A/" exp (V z ) |? 



(7) 



(iv) 



(5) 



where J\f is a normalization constant and taring) = 
X1/X2 for xi < X2- \ipo) is an entangled state of the 
spins for any value of 7 7^ 0. I.e. by varying 7 and 
choosing to, we have access to a rich variety of entangled 
many-particle states while staying in the lowest energy 
state. 

We proceed by discussing the conditions for adiabatic- 
ity and the sensitivity of the process to decoherence. In 
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this context particular attention has to be given to the 
phase transitions in cases (i) and (ii) since they are as- 
sociated with a merging of pairs of energies. Due to the 
symmetry of Q only one superposition of the associated 
states is however coupled to the nondegenerate initial 
ground state. Consequently for the transition to be adi- 
abatic it is sufficient that the characteristic time of the 
transfer T is much larger than the typical inverse fre- 
quency difference h/ AE to the next excited states. For 
the same reason collective decoherence processes caused 
e.g. by fluctuations in the external parameters £, A,Xi, 
and xi are suppressed in the present system by an expo- 
nential Boltzmann factor exp(— /3AE), where /3 _1 = ksT 
is the thermal energy of the heat bath. Decoherence pro- 
cesses caused by independent heat bath couplings of the 
individual spins do not have the symmetry of Q how- 
ever and hence do couple the degenerate states in cases 
(i) and (ii). Only in cases (Hi) and (iv), i.e. without a 
quantum phase transition, there is always a finite energy 
gap to all other states with the same total angular mo- 
mentum J. Thus the entanglement generation in cases 
(Hi) and (iv) will be robust against collective and indi- 
vidual decoherence processes, provided the energy gap is 
sufficiently large. 

It is not possible to give an analytic expression for the 
energy gap AE in the general case. Numerical investi- 
gations for up to 50 particles indicate that for a transfer 
efficiency close to 100% it is sufficient that 
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By choosing the trial function of the ground state close 
to the exact one it is possible to achieve 44 < 1, i.e. 
AEW ~ (-H")jv_2 - (H)n- For the simple trial func- 
tion I^W) = ai\m y = 0) + a 2 \m z = -N/2) one finds 
(H) N - 2 - (H) N = (3\ |£| XiX2|max, where (3 is a numer- 
ical factor of order unity, which varies only very slowly 
with N . Thus a reasonable estimate for the energy gap 
in case (Hi) is given by /3A |£| x|- This is also confirmed 
by our numerical calculations for particle numbers up to 
N = 50. If /3A |£| xi is sufficiently larger than the thermal 
energy of the environment, the probability of decoherence 
processes is strongly suppressed. 
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FIG. 2: Excitation scheme of pair of ions with ground state 
\g) and excited state \e) by bichromatic laser fields of equal 
and opposite detuning S. n denotes quantum number of trap 
oscillation. 




with A > N for case (i) and A > 1 for cases (ii — iv). 
T is the characteristic transfer time. An estimate for 
AE in cases (ii) and (Hi), i.e. for \i = can be ob- 
tained as follows: Using a variational method with a 
trial functions one finds an energy estimate (H)n 

for the ground state with J = N/2, E^ N) < (H) N . 
Secondly one can apply Temple's formula 0j to ob- 



tain a lower bound E. 



(N) 



> (H) 



N 



where 



(AH 2 )n = ((H — (H)) )n is the energy fluctuation in 

the trial state, e[ N ^ is the energy of the first excited 
state. One can show that Temple's formula gives the 



best lower bound when using only (AH 



(H)n and 



E) as parameters. Furthermore we make use of the in- . , , , ' , t i , a , n 

1 _ , N s (N-2) the two internal states and J z = 2 2^i=i( cr ee — ""go); wl ^h 



equality e[ N ^ > Eq N ~ 2 ^ between the energy of the first 
excited state for N particles and the ground state for 
N — 2 particles, in both cases with maximum J, which 
can easily be proven. Applying this inequality iteratively 
leads to 



Let us now discuss a possible implementation of the 
Lipkin Hamiltonian with fi — in an ion-trap sys- 
tem. Consider a linear trap with a string of ions with 
two relevant internal levels \g) and \e). The ions are as- 
sumed to be cooled such that only the in-phase collective 
oscillation of all ions is excited. The corresponding oscil- 
lation frequency is denoted by v. The two internal levels 
of the ions are coupled by two laser fields with frequen- 
cies Wi and u)2, and slowly- varying Rabi- frequencies 
and r^2- Assuming that both fields couple all ions in the 
same way, we can describe the system by the Hamilto- 
nian H = Ho + Hi n t, where Hq = hue) c + hu eg J z , c and 
c> being the annihilation and creation operators of the 
trap oscillation. huj eg is the energy separation between 

2 2^ j=l V" ee " gg> 

o~\iu = \l JL )ii(t L \ being the projector to the internal state 
of the ith ion. The interaction Hamiltonian H mt is 
given in rotating wave approximation by 



AffW > 1 



(ii) 



VI -AA 



where 



A > 



(AH^) N 



(9) 



(10) 



where rj is the Lamb-Dicke parameter. We assume that 
the laser frequencies have equal and opposite detuning 5 
from resonance wi^ — u eg ±5. 5 is assumed to be large 
compared to the linewidth of the resonance but suffi- 
ciently different from the frequency of the trap oscilla- 
tion, i.e. \5\, \ v±8\ 3> 7. As a consequence the dominant 



4 



processes are two-photon transitions leading to a simul- 
taneous excitation of pairs of ions as indicated in Figure 
121 For 0,% = Q2 this scheme has first been considered by 
S0rensen and Maimer 0, 0] in the context of quantum 
computation and dynamical entanglement generation. 
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FIG. 3: Numerical simulation of the adiabatic transfer for a 
system of four ions. Shown are populations in states \m y = 
0} and \m y = ±J) for fi 1]2 (i) = (a/2) (tanh(t/Ti, 2 ) + 1) as 
shown in the inset and for a/v = 0.6, vT\ = 2000, vT 2 = 1500, 
and 8/u — 0.9 (a), case (i); and S/u = 1.1 (b), case (iii) . 

We now assume that the ion trap is in the Lamb-Dicke 
limit, i.e. that the ions are cooled sufficiently enough, 
such that for all relevant excitation numbers n of the 
trap oscillation (n + l)rf <C 1 holds. In this limit one 
can expand the exponent in (|11|) to first order in 77. For 
large values of \S\ it is convenient to consider this inter- 
action in terms of a coarse-grained Hamiltonian which 
neglects the effects of rapidly oscillating terms. Using 
the time-averaging method of ref. one arrives at an 
effective Lipkin Hamiltonian JQ| with the identifications 
£ = (2iV )/(<5 2 - v 2 ), A = and X i,2 = «i T fl 2 . 

ill = ^2 corresponds to the case xi = M = m O which 
is exactly solvable and has been discussed in Refs.[^ll6|. 
It has also been shown in 0, ^| that the effective Hamil- 



tonian describes correctly the dynamics of the ions in 
a coarse-grained time picture. Furthermore the coupling 
scheme has successfully been implemented in experiments 
to generate entanglement between 4 ions Q. 

In Fig we have shown as an example the effective 
dynamics of a system of 4 ions driven by fields Qi and 
^2 for the cases (i) and (iii) following cq.Q) with /i = 0. 
One recognizes that a nearly perfect transfer is possible. 

In summary we have shown that it is possible to gen- 
erate specific entangled many-particle states in an en- 
semble of spins interacting through a collective coupling 
of the Lipkin-Meshkov-Glick type by adiabatic ground- 
state transitions. Two scenarios (i,ii) involve a two-fold 
degenerate ground state at some stages while two oth- 
ers (iii) and (iv) always have a nondegenerate ground 
state. In all cases there is a finite energy gap to other 
excited states. This gap can be rather large and thus fast 
processes are possible despite the required adiabaticity. 
In the asymptotic limits of the adiabatic transfer, the 
spectrum and eigenstates of the LMG-Hamiltonian can 
be exactly calculated. Furthermore for A = 1 there is a 
supersymmetry allowing for an explicit construction of 
the ground state for all times. Due to adiabaticity, all 
transitions are robust against parameter variations. Fur- 
thermore due to the symmetry of the coupling and the 
finite energy gap from the (degenerate or nondegenerate) 
ground state to excited states, collective decoherence pro- 
cesses are suppressed. In addition in the case of a non- 
degenerate ground state (case iii and iv) also individual 
decoherence processes are suppressed. The collective spin 
Hamiltonian can be implemented in a cold ensemble of 
ions in a linear trap driven by a nearly-resonant bichro- 
matic fields. 
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